Ellipsometry with an undetermined polarization state 
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We show that, under the right conditions, one can make highly accurate polarization-based mea- 
surements without knowing the absolute polarization state of the probing light field. It is shown 
that light, passed through a randomly varying birefringent material has a well-defined orbit on the 
Poincare sphere, which we term a generalized polarization state, that is preserved. Changes to the 
generalized polarization state can then be used in place of the absolute polarization states that make 
up the generalized state, to measure the change in polarization due to a sample under investigation. 
We illustrate the usefulness of this analysis approach by demonstrating fiber-based ellipsometry, 
where the polarization state of the probe light is unknown, and, yet, the ellipsometric angles of the 
investigated sample and A) are obtained with an accuracy comparable to that of conventional 
ellipsometry instruments by measuring changes to the generalized polarization state. 



I. INTRODUCTION 



Polarization-based measurements are increasingly im- 
portant, both as a fundamental tool for scientific re- 
search [1-8 , and as a vital tool for applications in the 
chemical, food, and pharmaceutical industries. Indeed, 
proposed classical communications systems, fluorescence 
measurements, non-invasive blood-sugar measurements, 
LIDAR, imaging through scattering media, strain, and 
temperature sensors all have implementations that rely 
on polarization based measurements j9HT3]. In surface 
science, it is well known that polarization based mea- 
surements are very sensitive. For instance, ellipsometry 
has been used to measure sub-monolayer changes in sur- 
face coverage [T4j|15]. As such, it is a nearly-ideal tech- 
nique for monitoring modern epitaxial fabrication tech- 
niques [16] and for contamination monitoring of vacuum 
components (e.g., synchrotron optics) [T7l|T8]. The sen- 
sitivity of ellipsometry comes from the ability to set and 
measure polarization states with a high degree of accu- 
racy. As a result, it is generally thought that, in order to 
make polarization-based measurements, one must know 
the polarization state of the probing light field through- 
out the optical train, and, especially, the polarization 
state of the light incident on the surface of interest. The 
corollary to this is that the use of optical components 
that disturb the polarization (e.g., optical fibers) require 
calibration, and, in the case of optical fibers, where tem- 
perature and stress induced birefringence have a large 
influence on the polarization state, this is not possible 
in all but the most limited circumstances. This is be- 
cause, in contrast to what the name suggests, polariza- 
tion maintaining (PM) fibers do not preserve arbitrary 
polarization states [19^. For instance, if linearly polar- 
ized light is injected at an arbitrary polarization orienta- 
tion with respect to the axes of birefringence of the fiber, 
then the output is an elliptically polarized state with an 
unpredictable orientation and degree of ellipticity. As 
a consequence, fiber-based ellipsometry and polarimetry 
have relied on rather complicated experimental appara- 



tus [201 EI] 7 use a wavelength that is only supported by 
one polarization mode in a PM fiber [22 , or only provide 
qualitative information [23j . 



The core discovery that we present in this paper is 
that, given an input polarization state of light and an 
optical fiber that is subject to environmentally induced 
birefringence variations, the output state does not map 
to every point on the Poincare sphere, but rather, only 
a discrete set of polarization states are accessible. As 
the temperature, for instance, varies in time, the set of 
output polarization states appear as a fixed orbit on the 
surface of the Poincare sphere. In our work, we show that 
the existence of such orbits, as an intrinsic property of the 
fiber, can be seen as a higher dimensional type of polar- 
ization preservation that survives severe environmental 
perturbances. Once this orbit is known, polarization- 
based measurements can be made by measuring devia- 
tions from the orbit. This is very different from direct 
calibration, where an input polarization state is mapped 
to a single output polarization state for a known set of 
environmental parameters. Instead, a single input po- 
larization state is mapped to a set of output states as 
the environment varies over some (generally unknown) 
range of temperatures and stresses. This approach turns 
a measurement problem — environmental noise — into an 
advantage by using the statistical properties of the noise 
to provide increased sensitivity. 



In this paper, we present experimental data and mod- 
eling results that demonstrate and make use of the obser- 
vation that the polarization follows a random path, but, 
importantly, remains bound to well-defined subspace. As 
an example of how such knowledge can be applied, we 
demonstrate a fiber-based ellipsometer. We show that 
this ellipsometer is capable of detecting a well-defined 
carbon layer that is less than 1 nm thick on top of a 
multilayer Bragg reflecting mirror. 
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FIG. 1: Schematic of the eUipsometer setup. The polariza- 
tion control optics (PC) set the light's input polarization, as 
illustrated on the Poincare sphere (a). After passing through 
a non-polarizing beam splitter (BS), part of the light is sent 
to a rotating compensator for analysis (PM 1) of the input 
polarization and part is coupled into the PM fiber. The polar- 
ization state of the light exiting the fiber falls somewhere on 
the orbit illustrated on Poincare sphere (b) . The light reflects 
off the sample, rotating the orbit to that shown in Poincare 
sphere (c) and the final polarization state is measured at PM 
2. 



II. EXPERIMENTAL SETUP 

The experimental setup is shown schematically in 
Fig. [l] As a light source, we used a Helium-Neon laser 
that emitted light that was nominally linearly polarized. 
The polarization of the light field was purified to a lin- 
early polarized state, of which the azimuth angle was set 
to be parallel to the plane of the optical table, using a 
Wollaston prism. The Wollaston prism ensures that the 
polarization purity is 100000:1. The input polarization 
state to the fiber was controlled using a Soleil-Babinet 
compensator or a quarter- wave plate, which were man- 
ually controlled using precision rotation and translation 
stages. This allowed the individual relative Stokes vec- 
tor intensities to be set within 0.02% of a desired value. 
After setting the polarization, the laser beam was passed 
through a non-polarizing beam splitter, before being cou- 
pled into a 2 m long PM-fiber (Thorlabs, PM-630-HP) 
using a glass aspheric lens. Light was coupled out of the 
output end of the fiber and collimated with an identical 
coupling lens. The output fiber end was mounted on a 
manual rotation stage to set the angle of incidence with 
respect to the sample's plane of reflection. The accuracy 
of the angle of incidence on the sample was found to be 
0.5°. In all measurements, the middle section of the fiber 
(about 1.3 m) was coiled up on a metal cylinder and 
placed in a water bath with a heater to simulate envi- 
ronmental changes by varying the temperature between 
19 and 30° C. Only a small fraction of this range (2-3° C) 
was used during the experiments, though, because that 
was sufficient to obtain the full range of accessible polar- 
ization states (see below). 



After reflection from the sample, the polarization state 
of the light was measured using a rotating compensator, 
which consisted of a quarter waveplate, rotating at an 
angular frequency of 0.1 Hz, a Wollaston prism, and a Si 
PiN photodiode (PM 1 and 2 in Fig. [!}. The photodi- 
ode voltage was observed on an oscilloscope and trans- 
ferred to a computer for analysis. The photodiode signals 
were evaluated to obtain the polarization in terms of the 
Stokes vector components The polarization states 

are calculated as Stokes vectors and presented graphi- 
cally in either 2D or 3D projections of the Stokes vector 
space. In order to restrict our analysis to the polarization 
properties of the PM fiber and sample, we compensate for 
the (small) modification of the polarization state due to 
the non-polarizing splitter with a predetermined matrix. 

Ideally, neither the PM fiber nor the sample are chang- 
ing, and the frequency components of the photodiode sig- 
nal should only contain the DC component and two har- 
monics of the rotating compensator's angular frequency. 
However, this is true only when the polarization state is 
constant during the measurement period. In our case, 
the fiber's output states are constantly changing and, 
thus, allows a small non-zero amplitude for additional 
frequency components. By analyzing the amplitude of 
these components, the direct polarization state measure- 
ment noise was estimated to be 0.02° (for the ellipsomet- 
ric parameter A), which is comparable to a conventional 
EUipsometer [25] . 



III. RESULTS 

The measured polarization states of the output of the 
PM fiber are plotted in Fig. 2|a) as zero-dimensional 
points on the two-dimensional surface of the Poincare 
sphere, where all data are, unless it is specifically stated 
otherwise, normalized so that the total intensity is unity. 
The notable feature is that despite the wide variation 
in temperature, the polarization states are not observed 
to fall randomly on the two-dimensional Poincare sur- 
face. Instead, a single, one-dimensional orbit is traced 
out, which is shown more clearly by taking a cross sec- 
tion through the sphere to obtain a projection of the 
orbit in a Stokes vector plane (Fig. |2|b)). The analy- 
sis of the orientation of this orbit is key to performing 
polarization-based measurements without precise knowl- 
edge of the polarization of the probing light field. 

To show that this polarization orbit can be used as a 
generalized, higher dimensional polarization state and be 
used to make polarization-based measurements, we show 
how the two ellipsometric parameters, ^ and A can be 
calculated from changes to the polarization orbit. We 
begin by considering light exiting a fiber in a set of po- 
larization states denoted by X^^ where the subscript and 
superscript denote the polarization coordinate systems 
of the sample {s — p) and the fiber {x — y). An isotropic 
reflecting surface transforms the polarization state of in- 
cident light according to the following Mueller matrix in 
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FIG. 2: Polarization state of the light exiting the fiber for 
a fixed input polarization and varying the fiber temperature 
between 19.4 and 21.8°C. Subfigure (a) shows a three dimen- 
sional image of the polarization state on the Poincare sphere, 
while (b) is the projection of that orbit onto the 3rd and 4th 
Stokes vector components. The solid line is a fit to the data 
with an ellipse (see below). 



the sample coordinate system [25] 



M = A- 
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^ /)/2 and Ir^^l and Ir^yl are inten- 
sity reflectivities for the P and S polarized light respec- 
tively. The initial states become 



where A 



(2) 



In other words, the ellipsometric information of the sam- 
ple is carried in a global transformation from XJ^ to IZ^^ . 

The matrix M is a direct sum of the matrices of its 
two block diagonal subspaces. The (desired) values of ^ 
and A can, thus, be obtained separately. 



The set of polarization states of the light field incident 
on, and reflected by the sample are given by 
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substituting equations |3] into equation [2] gives 

Ri = A{Ii - h cos 2^) 

R2 = A{l2 - h cos 2^) 

Rs = Asin(2^)(/3Cos A - /4sin A) 

R4 = Asin(2^)(/3sinA + /4CosA) 

Dividing Eqs. |4] by [5] and rearranging gives: 



^ = - arccos 
2 



1 



Rr 



(4) 
(5) 
(6) 
(7) 



(8) 



where Rr = R2/R1 and Ir = l2/h- A is obtained by 
requiring consistency between equations [6] and ^ Note 
that equations |6] and [7| have the form of a scaling factor 
and a rotation, which, when applied to an ellipse, reduces 
its area and changes the orientation of its major axis in 
Stokes space. In practice, since the polarization orbit, 
when projected in the plane of the 3rd and 4th Stokes 
vectors, is an ellipse, A is the angle difference between 
the orientations of the ellipses given by XJ^ and TZ^^ . 

It is important to note that our analysis assumes that 
all states are measured in the x — y coordinate system. 
Slight changes of alignment between the coordinate sys- 
tem of the fiber and the sample introduces an extra ro- 
tation to the polarization state that depends on the mis- 
alignment angle, 5a. In that case, it can be shown that 
the systematic error added to ^ and A is of the order 
of Sa'^. Provided that the physical alignment of the PM- 
fiber is held sufficiently constant, the uncertainty in A 
and ^ should be comparable to that of conventional el- 
lipsometry. 

To demonstrate the applicability of this analysis, we 
performed fiber-based ellipsometry on three multilayer 
Bragg reflecting samples [26 , two of which have been 
coated with amorphous hydrogenated carbon [171 HE]- 
Spectroscopic ellipsometry (Wollam M2000) revealed 
that the carbon layers on these samples are 0.3, and 
0.8 nm thick, while the third, supposedly uncoated sam- 
ple, has ~0.1 nm naturally occurring carbon layer with 
a different composition to that of the first two samples. 
Both conventional and fiber-based ellipsometry measure- 
ments were performed at an angle of incidence of 66°. 
The raw ellipsometric data from the fiber-based ellip- 
someter, using the experimental procedure described de- 
scribed above, is shown in Fig. [sja). The measured or- 
bits have been projected onto the plane of the third and 
fourth Stokes vector components (the inset shows the 
three dimensional representation of the data). As illus- 
trated by the zoomed in section, shown in Fig. |3|b), the 
data for each sample are systematically modified by the 
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FIG. 3: Subfigure (a): projection of the polarization orbits on 
the plane of the third and fourth Stokes vector components for 
the light exiting the fiber (solid green circles) , after reflection 
from: a multilayer mirror (filled black diamonds), a multilayer 
mirror with a 0.3 nm thick carbon layer (open red circles), and 
a multilayer mirror with a 0.8 nm thick carbon layer (open 
blue diamonds). Subfigure (b) shows a zoomed section of the 
projection of the polarization orbits. The lines are the fitted 
ellipses, which are used to determine ^ and A. 

sub-nanometer layers of carbon to generate three differ- 
ent ellipses. The data for each sample show very little 
deviation from the fitted ellipse. 

Using the analysis described above, we obtain the ^ 
and A values presented in Tab. [l[ These values agree 
well with those obtained from a conventional ellipsome- 
ter. The small systematic difference between the values 
from the two instruments is due to a small difference 
in the angle of incidence between the two measurements 
(the accuracy of the angle of incidence is ±0.5°). 

IV. DISCUSSION 

To understand why PM fibers that are subject to vary- 
ing temperature and strain, produce a polarization orbit. 



TABLE I: ^ and A values for fiber-based ellipsometric mea- 
surements (columns 2 and 3) and conventional ellipsometer 
measurements (columns 4 and 5) on carbon-coated MLMs 



Sample 


^ ° 


A ° 




A ° 


Bare MLM 


28.04 


152.5 


28.158 


152.77 


0.3 nm 


28.00 


152.1 


28.265 


152.39 


0.8 nm 


28.48 


148.9 


28.389 


150.81 



a numerical model of a non-ideal PM fiber was devel- 
oped. In simple terms, the fiber is a strongly birefringent 
material with a time and space-dependent beat length, 
and some degree of polarization mode cross-talk. In this 
model, the orientation of the axes of birefringence and the 
refractive index difference between the two axes changes 
as a function of position along the fiber and time. To an- 
alyze the effect of such a fiber, we break the fiber up into 
ideal segments, Al. Mathematically, this corresponds to 
rotating the coordinate system of the fiber for light en- 
tering a fiber segment and then reversing the rotation 
upon exiting. The phase delay experienced by the light 
field then depends on its orientation relative to the ro- 
tated coordinate system. This introduces a polarization 
change that is purely due to the relative phase delay of 
the components of the electric field along the local axes 
of birefringence. The transfer matrix of the entire fiber 
is then given by a series of matrix multiplications [25] : 

Mf, = n^(Aa(/,t))m^(^)i?(-Aa(/,t)) (9) 

R{Aa{l^t)) and R{—Aa{l^t)) are rotations due to an an- 
gular (Aa) misalignment of the axes of birefringence of 
two neighboring ideal segments, while is the matrix 
transfer function due to the phase delay. The phase delay 
is expressed as a function of the beat length, 6(/,t), nor- 
malized to the segment length. Both b{l^t) and Aa{l^t) 
are time dependent in a constrained but random man- 
ner. In these simulations, the maximum limits of Aa 
and changes to b{l^t) were given by the yield strength of 
glass — the strain required to introduce these variations 
should not break the fiber. To allow direct comparison 
between the model results and experimental results, the 
average variation in Aa was determined by specifying a 
maximum polarization cross-talk per unit length. 

As the phase change varies as a function of position 
and time in the fiber, the output polarization states form 
a simple, circular orbit in Stokes space. The effect of 
cross-talk is to introduce deviations from a simple orbit 
in Stokes space, instead, a more complicated, but repeat- 
ing, orbit is obtained. Numerically evaluating equation 
|9] under conditions that the maximum total cross talk 
is 20 dB, which is the maximum specified cross-talk for 
our fiber, shows (see Fig.[4| that, even in the presence of 
a small amount of cross-talk, the polarization orbit can 
still be approximated by a simple orbit (the elongation 
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FIG. 4: Numerical evaluation of equation [9] for a total polar- 
ization mode cross talk of 20 dB. Time-dependent tempera- 
ture changes cause changes to the local degree and orientation 
of the fiber's birefringence. The resulting output is an orbit 
on the Poincare sphere. 



avoided by making measurements at more than one angle 
of incidence. 

In conclusion, we have demonstrated that sensitive 
polarization-based measurements are possible even when 
the polarization state of the probing light field has been 
passed through an optical element with an unpredictably 
fluctuating birefringence. This is possible because, for 
a fixed input polarization, the environmentally induced 
fluctuations in the output polarization after propaga- 
tion through an optical element or material still lie in 
a single orbit of the Poincare sphere, creating what we 
have termed a generalized, higher dimensional polariza- 
tion state. Polarization-based measurements can then 
be made by analyzing changes to the higher dimensional 
state, rather than the individual polarization states that 
make up the higher dimensional state. To illustrate this, 
we demonstrated a fiber-based ellipsometer, capable of 
detecting carbon layers with a thickness of 0.3 nm. 



of the red markers in Fig. [4] indicate the deviation from 
a simple orbit). Importantly, these simulations indicate 
that, as long as we can quantify the polarization orbit of 
an optical component, then it will be possible to use light 
exiting that component for polarization-based measure- 
ments. It should be noted there is the (slim) possibility 
that the polarization change induced by the sample is a 
rotation in the plane of the polarization orbit. In this 
case, the sample is undetectable. This problem can be 
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